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Abstract. We report on the development of a new model for the underlying event in hadron-hadron 
collisions. The model includes parton showers for all interactions, as well as non-trivial flavour, momentum, 
and colour correlations between interaction initiators and beam remnant partons. 
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1 Introduction 

A simple consequence of the composite nature of hadrons 
is the possibility to have hadron-hadron collisions in which 
several distinct pairs of partons collide with each other 
(multiple interactions). In fact, simple perturbative calcu- 
lations can be used to show [1] that most inelastic events in 
high -energy hadronic collisions should contain several per- 
turbatively calculable interactions, in addition to what- 
ever nonperturbative phenomena may be present. 

Although most of this activity is not hard enough to 
play a significant role in the description of high-j>j_ jet 
physics, it can be responsible for a large fraction of the 
total multiplicity (and large fluctuations in it), for semi- 
hard (mini-)jets in the event, for the details of jet profiles 
and for the jet pedestal effect, leading to random as well 
as systematic shifts in the jet energy scale. Thus, a good 
understanding of multiple interactions would seem pre- 
requisite to carrying out precision studies involving jets 
and/or the underlying event in hadronic collisions. 

In an earlier study [1], it was argued that all the un- 
derlying event activity could be explained by the multiple 
interactions mechanism alone. However, while the origin 
of underlying events is thus assumed to be perturbative, 
many nonperturbative aspects still force their entrance on 
the stage. This in particular relates to the structure of 
beam remnants and to the correlations in flavour, colour, 
and momentum between the partons involved. In [1], only 
very simple beam remnant structures could technically be 
dealt with, hence substantial simplifications had to be im- 
posed. 

In recent years, the physics of the underlying event 
has come to attract more attention. Simple parameteri- 
zations can be tuned to describe the average underlying 
activity, but are inadequate to fully describe correlations 
and fluctuations. The increased interest and the new data 



now prompts us to develop a more realistic framework for 
multiple interactions than the one in ref. [1] , while making 
use of many of the same underlying ideas. 

One new aspect was the augmentation in [2] of the 
standard Lund string fragmentation framework [3] to in- 
clude the hadronization of colour topologies containing 
non-zero baryon number. In the context of multiple in- 
teractions, this improvement means that almost arbitrar- 
ily complicated baryon beam remnants may now be dealt 
with, hence many of the restrictions present in the old 
model are no longer necessary. 

Here, we present a model for how flavours, colours, and 
momenta are correlated between all the partons involved 
in a hadron-hadron collision, both those that undergo 
interactions and those constituting the beam remnants. 
However, all aspects of the model cannot be treated within 
the limits of this format, hence some aspects have been 
left out; we concentrate exclusively on baryon beams and 
neither the dependence on impact parameter nor the as- 
signment of 'primordial k±' to parton shower initiators is 
addressed here. More complete descriptions may be found 
in [4,5]. 

This article is organized as follows. In Section 2 the 
main work on flavour and momentum space correlations 
is presented and in 3 the very thorny issue of colour cor- 
relations. Finally, Section 4 provides a brief summary. 



2 Towards a Realistic Model 

Consider a hadron undergoing multiple interactions in a 
collision. Such an object should be described by multi- 
parton densities, giving the joint probability of simulta- 
neously finding n partons with flavours /i,...,/„, car- 
rying momentum fractions Xi,...,x n inside the hadron, 
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when probed by interactions at scales Q 2 , . . . , Q\. How- 
ever, we are nowhere near having sufficient experimental 
information to pin down such distributions. Therefore, and 
wishing to make maximal use of the information that we 
do have, namely the standard one-parton-inclusive parton 
densities, we propose the following strategy. 

As described in [1], the interactions may be generated 
in an ordered sequence of falling pj_ . For the hardest inter- 
action, all smaller p±^ scales may be effectively integrated 
out of the (unknown) fully correlated distributions, leav- 
ing an object described by the standard one-parton dis- 
tributions, by definition. For the second and subsequent 
interactions, again all lower-pj_ scales can be integrated 
out, but the correlations with the first cannot, and so on. 

Thus, we introduce modified parton densities, that cor- 
relate the i'th interaction and its shower evolution to what 
happened in the i — 1 previous ones. 

The first and most trivial observation is that each 
interaction i removes a momentum fraction Xi from the 
hadron remnant. Already in [1] this momentum loss was 
taken into account by assuming a simple scaling ansatz 
for the parton distributions, f(x) — ► f(x/X)/X, where 
A = 1 — Y^i=i x i i s the momentum remaining in the beam 
hadron after the n first interactions. Effectively, the PDF's 
are simply 'squeezed' into the range x G [0, A]. 

Next, for a given baryon, the valence distribution of 
flavour / after n interactions, qf vn (x,Q 2 ), should follow 
the counting rule: 
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qfvn{x,Q 2 ) Ax = N fvn , 



(1) 



where Nf vn is the number of valence quarks of flavour 
/ remaining in the hadron remnant. This rule may be 
enforced by scaling the original distribution down, by the 
ratio of remaining to original valence quarks Nf vn /Nf v o, 
in addition to the x scaling mentioned above. 

Also, when a sea quark is knocked out of a hadron, it 
must leave behind a corresponding antisea parton in the 
beam remnant. We call this a companion quark. In the 
perturbative approximation the sea quark q s and its com- 
panion q c come from a gluon branching g — > q s + q c (it is 
implicit that if q s is a quark, q c is its antiquark). Start- 
ing from this perturbative ansatz, and neglecting other 
interactions and any subsequent perturbative evolution of 
the q c distribution, we obtain the q c distribution from the 
probability that a sea quark q s , carrying a momentum 
fraction x s , is produced by the branching of a gluon with 
momentum fraction y, so that the companion has a mo- 
mentum fraction x = y — x s , 



q c {x, x s 



= C [ g(y)E 
Jo 

= C g(y) Pg^q s q. 
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with P B 



♦q B q< 



the usual DGLAP gluon splitting kernel and 



C a normalization constant which can be obtained by im- 



posing the counting rule: 
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q c (x;x s ) dx = 1. 



(3) 



The exact form of C depends on the shape assumed for 
the gluon disitribution. Qualitatively, however, any falling 
gluon distribution oc 1/x convoluted with the almost flat 
g — > qq splitting kernel yields companion distributions 
which tend to a constant q c (x;x s ) ~ C/2x 2 below x s 
and which exhibit power-like fall-offs q c (x;x s ) oc 1/x 2 
above it, with some modulation of the latter depending 
on the large-x behaviour assumed for g(x). Also note that 
xq c (x;x s ) should be peaked around x « x s , by virtue of 
the symmetric P g ^ qsqc splitting kernel. 

Without any further change, the reduction of the va- 
lence distributions and the introduction of companion dis- 
tributions, in the manner described above, would result in 
a violation of the total momentum sum rule, 
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Qfn( X > Q 2 ) + 9n(x, Q 2 ) ) dx 



A, 



(4) 



since by removing a valence quark from the parton dis- 
tributions we also remove a total amount of momentum 
corresponding to (xf v ), the average momentum fraction 
carried by a valence quark of flavour /: 



J X xq fvn (x,Q 2 ) dx 



(Xfy 



J X qf vn {x, Q 2 ) dx 



= X (x fvQ ) , (5) 



and by adding a companion distribution we add an anal- 
ogously defined momentum fraction. 

To ensure that eq. (4) is still respected, we assume 
that the sea+gluon normalizations fluctuate up when a 
valence distribution is reduced and down when a compan- 
ion distribution is added. In addition, the requirement of 
a physical x range is of course still maintained by 'squeez- 
ing' all distributions into the interval x E [0, A]. 

The full parton distributions after n interactions thus 
take the forms: 



Qfn {x) = 



N 



fvn 



Nfvo 
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+ 
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gn{x) = -|<?0 (^) > 



(6) 
(7) 



where we have suppressed the dependence on Q 2 for brevity, 
9/vO (<Z/so) denotes the original valence (sea) distribution 
of flavour /, and the index j on the companion distribu- 
tions qf Cj counts different companion quarks of the same 
flavour, /. As already mentioned, the normalization fac- 
tor a multiplying the gluon and sea distributions can be 
determined from overall momentum conservation in the 
incoming hadron: 



a = 



1 ~ Y,f N fvn(x fv o) - J2f,j(xf Cj o) 

1 -T,f N fv0(Xf v0 ) 



(8) 
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# Parton in hadron remnant 




Fig. 1. a) The initial state of a baryon, consisting of 3 va- 
lence quarks conneted antisymmetrically in colour via a cen- 
tral 'string junction', J. b) Example of how a given set of par- 
ton shower intitators could have been radiated off the initial 
configuration, in the case of the 'purely random' correlations 
discussed in the text. 



After the perturbative interactions have taken each 
their share of longitudinal momentum, the question arises 
how the remaining momentum is shared between the beam 
remnant partons. Here, valence quarks receive an x picked 
at random according to a small-Q 2 valence-like parton 
density, while sea quarks must be companions of one of 
the initiator quarks, and hence should have an x picked 
according to the q c (x;x s ) distribution introduced above. 
In the rare case that no valence quarks remain and no sea 
quarks need be added for flavour conservation, the beam 
remnant is represented by a gluon, carrying all of the beam 
remnant longitudinal momentum. 

Further aspects of the model include the possible for- 
mation of composite objects in the beam remnants (e.g. di- 
quarks) and the addition of non-zero primordial k± val- 
ues to the parton shower initiators. Especially the latter 
introduces some complications, to obtain consistent kine- 
matics. Details on these aspects will be presented in [5]. 



3 Colour Correlations 

The initial state of a baryon may be represented by three 
valence quarks, connected antisymmetrically in colour via 
a central junction, which acts as a switchyard for the 
colour flow and carries the net baryon number. This situ- 
ation is illustrated in Fig. la. 

The precise colour-space evolution of this state into the 
initiator and beam remnant partons actually found in a 
given event is not predicted by perturbation theory, but is 
crucial in determining how the system hadronizes; in the 
Lund string model [3], two colour-connected final state 
partons together define a string piece, which hadronizes 
by successive non-perturbative breakups along the string. 
Thus, the colour flow of an event determines the topology 
of the hadronizing strings, and consequentially where and 
how many hadrons will be produced. 

For the perturbative parts of the event, a unique colour 
flow may be consistently assigned [6], but for the connec- 
tions among initiator and beam remnant partons, addi- 
tional assumptions are necessary. The question can essen- 
tially be reduced to one of choosing a fictitious sequence 



of gluon emissions off the initial valence topology, since 
sea quarks together with their companion partners are as- 
sociated with parent gluons, by definition. 

The simplest solution is to assume that gluons are at- 
tached to the initial quark lines in a random order, see 
Fig. lb. If so, the junction would rarely be colour con- 
nected directly to two valence quarks in the beam rem- 
nant. It should be clear that the migration of the baryon 
number depends sensitively upon which partons in the fi- 
nal state the junction ends up being connected to (see [2] 
for details on junction fragmentation). Thus, if the connec- 
tions are purely random, the baryon number of the initial 
state should quite often be disconnected from the beam 
remnant altogether and be able to migrate to both large 
p± and small xf values. Empirically, this may not be de- 
sireable, hence a free suppression parameter is introduced 
to suppress gluon attachments onto colour lines that lie 
entirely within the remnant. 

Finally, we imagine a few different possibilities for the 
ordering of the emissions off a given colour line: 1) ran- 
dom, 2) gluons are ordered according to the rapidity of 
the hard scattering subsystem they are associated with 
(where beam remnant partons are assigned a fixed large 
rapidity in the direction of their parent hadrons), and 3) 
gluons are ordered so as to give rise to the smallest possible 
total string lengths in the final state. The two latter possi- 
bilities represent different attempts to minimize the total 
potential energy of the system (the string length), since 
the former seems to result in a too large hadron multiplic- 
ity per interaction. A variable which we have found to be 
sensitive to the choice of colour topology is the average p± 
vs. n c hargcd, but our studies are not concluded yet. Thus, 
we do not pretend to have the final solution to these ques- 
tions. Rather, the colour correlations, both in the initial 
and in the final state, still represents the major open issue 
in our studies. 

4 Conclusion 

The development of a new model for the underlying event 
in hadron-hadron collisions has been reported. This model 
extends the multiple interactions mechanism proposed in 
[1] with the possibility of non-trivial flavour and momen- 
tum correlations, parton showers for all initiator and final 
state partons, and several options for colour correlations 
between initiator and beam remnant partons. Many of 
these improvements rely on the development of junction 
fragmentation in [2] . 

The issue of colour correlations is still actively under 
study. 
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